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FOREWORD 


This report covers the work completed on the research project "Analysis 
and Computation of Internal Flow Field in a Scramjet Engine." The work was 
supported by the NASA Langley Research Center (Computational Methods Branch 
of the High-Speed Aerodynamics Division) through research grant NAG-1 -423, 
and monitored by Dr. Ajay Kumar of the High-Speed Aerodynamics Division. 
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A SECOND-ORDER ACCURATE PARABOLIZED 
NAVIER-STOKES ALGORITHM FOR INTERNAL FLOWS 


By 

Tawit Chitsomboon 1 and S. N. Tiwari 2 
SUMMARY 

A parabolized Navier-Stokes algorithm is presented which is implicit 
and of second-order accuracy in both the cross flow and marching directions. 
The algorithm is used to analyze three model supersonic flow problems (the 
flow over a flat plate with an impinging shock, and inviscid and viscous 
flows over a 10-degree wedge). The results are found to be in good agree- 
ment with the results of other techniques available in the literature. 

1. INTRODUCTION 

Parabolized Navier-Stokes (PNS) equations have recently become popular 
in solving steady supersonic flows. The increase in the use of the PNS 
equations is due to the fact that the algorithm is very efficient with 
regard to both execution times and storage requirements. 

Unlike boundary-1 ayer type equations, the PNS equations allow cross- 
stream interaction because the normal momentum equation is retained. Rever- 
sal of the flow field in the streamwise direction, however, is not permitted 
because the parabolic nature of the equations dictates no upstream influ- 
ence. 

In order for the code described in this report to be applicable, the 
flow field must be supersonic. A small fraction of a subsonic layer close 
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to solid walls is permitted, but the layer must not be too big as to domi- 
nate the entire flow field. These conditions usually are satisfied for high 
Reynolds number supersonic flows over bodies of small geometrical variations 
in the streamwise direction. 

The algorithm presented in this study is implicit and of second-order 
accuracy in both cross flow and marching directions. The generalized trans- 
formation of coordinates is used for convenience of boundary-conditions 
implementation. The algorithm is used to analyze a model two-dimensional, 
supersonic inlet problem at high Mach number and a supersonic flow over a 
flat plate with an incident shock wave. The results of these model problems 
are compared with the results of other techniques. 

2. THE PNS EQUATIONS FOR TWO-DIMENSIONAL FLOW 

The governing equations for the present study are presented here brief- 
ly. For a descriptive representation of both two- and three-dimensional 
equations, reference should be made to Refs. 1 and 2. 

The Navier-Stokes equations for a steady two-dimensional flow without 
body forces and internal heat generation can be written in a non-dimensional 
conservation-law form (for a Cartesian coordinates system) as 

— (E - E ) +_L (F - F ) - 0 (2.1) 
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To close the system of equations, the following perfect gas equation of 
state is used. 


p ■ (Y - 1) pe (2.2) 

The coefficient of viscosity is determined from Sutherland's model. 

Next, equation (2.1) is transformed into the body-conforming coordi 
nates and put into a strong conservation-law form as follows: 

_i (E) + JL ( F) = 0 (2.3) 

9S 9n 


where 

E = -i [5 X E + 5 y F] 

F = — [n x (E - E v ) +n y (F - F v )] 

J 

and J represents the Jacobian of transformation, and 5 V , n v , n.. 

a y a y 

are the metrics of transformation. In order to put equation (2.3) in the 
PNS form all the viscous cross-derivative terms in F have been neglected. 

Unless properly treated, the streamwise pressure gradient in the sub- 
sonic layer will cause the marching procedure to be divergent. In this 
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study, Vigneron's technique, in which only a fraction of the streamwise 
pressure gradient is retained (Ref. 2), is used as opposed to the sublayer 
approximation employed in Refs. 1 and 5. 

The modified equation suitable for a stable marching is then expressed 
as: 


ill + il_ + if =o (2.4) 

35 3? 3n 


P = — [0, (l-w)'5 x P, (l-w)5yP> 0] T 
J 

where w is a parameter obtained from an eigenvalue analysis (Ref. 2). 

3. FINITE-DIFFERENCE ALGORITHM 

A fully implicit, nori-iterati ve, finite-difference algorithm, as de- 
scribed by Schniff and Steger (Ref. 1), is employed in this study. The 
resulting algorithm is in delta form and requires a block-tridiagonal in- 
version at each marching step. Hie details of the algorithm are available 
in Ref. 1 and the final form is expressed here as: 


(3.1) 


[A 1 + (1 - a) A 5 6 B 1 ] A n q = 

- (V - A 1 " 1 ) q 1 * + a (E* 1 - t Ul ) 

- (1 - a) A?6 n F 1 - A 1 -1 P 
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In equation (3.1) If a a 0 the algorithm is first-order accurate and 
if a 3 1/3 It Is second-order accurate. The quantities A and 8 are 
the Jacobian matrices resulting from the linearization of E and F with 
respect to q, respectively; they are 4x4 matrices for the two-dimen- 
sional case. The tildes above A, B and F indicate- that the metric 

coefficients are to be evaluated at the station (i+1), whereas the Jacobian 

h 

of transformation is evaluated at the i station. 

For the viscous case, the boundary conditions are two no-slip condi- 
tions, zero normal pressure gradient, and a specified temperature at the 
solid wall. The method of simple extrapolations has been used successfully 
to implement the boundary conditions for the inviscid cases. For this 
method to work well, it is necessary to cluster the grid points near the 
solid wall. 

With an appropriate initial data plane, the PNS algorithm of equation 
(3.1) can be marched in the ? direction in one sweep to obtain a solution. 
Because this is an implicit- scheme, the solution procedure requires an in- 
version of a block-tridiagonal system. 

4. NUMERICAL SMOOTHING 

For viscous as well as inviscid marching, both implicit and explicit 
smoothing terms of second and fourth order, respectively, are employed to 
dampen the oscillation. The implicit smoothing function also serves to make 
the diagonal terms of the block-tridiagonal system non-zero. 

By numerical experiments, it is found that smoothing functions on q 
vector are more effective than smoothing functions on the flux vector E 
which is used in Refs. 1 and 5. 

The second-order implicit smoothing term of the form 


is added to the left-hand side of equation (3.1), whereas the fourth-order 
explicit smoothing term is implemented as 

- -1 (VA)2 J 1 q 1 

and is added to the right-hand side of equation (3.1). 

From linear stability analysis, the upper-bound value of e £ is 1/8. 
For the present model problems, the explicit smoothing function alone cannot 
stabilize the marching procedure. The marching will be stable if the value 
of £j is bigger than a value which is proportional to the step size. This 
value of ej has to be determined numerical ly. It is found also that if 
the value of e T is two to three times that of e_, the upper-bound value 
of e E is not restricted. 


5. RESULTS AND DISCUSSION 

Three test cases have been run to test the PNS algorithm. A brief dis- 
cussion of the results is presented here. 

5.1 Inviscid Flow with an Incident Shock Wave on a Flats Plate 

The schematic diagram for this case is illustrated in Figure 5.1. The 
initial data plane was obtained from the exact solution of a Mach-2.9 flow 
over an 11-degree wedge. 

A pressure-contour plot of the solution is shown in Figure 5.2. Figure 
5.3 compares the normalized surface pressure on the flate plate with the 





Exact Solution 
Inviscid Marching 



Figure 5.3. Surface pressure distribution for supersonic flow over a flate plate. 


exact solution obtained from a shock table. Hie pressure at 1/4 inch from 
the plate is compared with the MacCormack's predictor-corrector method and 
Osher's flux-splitting method in Figure 5.4. The results are seen to be in 
very good agreement with the results of other techniques. 

5.2 Inviscid Internal Flow Over a 10-Degree Wedge, = 5 

The geometry and free stream conditions of this case are shown in Fig- 
ure 5.5. The initial data plane for this case is specified as a uniform 
flow of free-stream conditions at 1 cm. from the wedge. 

Figure 5.6 shows a pressure-contour plot of the solution. The results 
for surface-pressure distribution shown in Figure 5.7 indicate a very good 
agreement with the results of MacCormack's method of Ref. 3. The pressure 
distribution at 2 cm. from the solid wall is compared with that of Ref. 3, 
using the MacCormack's method, in Figure 5.8. It is seen that the peak 
pressure of the inviscid-marching solution is delayed downstream by almost 
1 cm. More investigations are needed to improve the agreement. The oscil- 
lation after the peak pressure for the MacCormack solution is non-physical 
and is due to the grid-size aspect ratio. 

5.3 Viscous Internal Flow Over a 10-Degree Wedge, » 5 

The geometry and free-stream conditions of this case are the same as 
for the inviscid case. The initial data plane is specified at 1 cm. from 
the wedge and is ootained from the solution of Ref. 3. 

The pressure-contour plot in Figure 5.9 shows a splitting of the shock 
wave; this is due to the oscillation across the shock wave. Figure 5.10 
indicates a fairly good agreement in the surface pressure with the full 
Navier-Stokes solution of Ref. 3. The pressure at 2 cm. from the wall is 
compared with that of Ref. 3 in Figure 5.11. As in the inviscid case, the 
peak pressure is seen to be delayed in the downstream region. 
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Figure 5.4. Pressure distribution at 1/4 inch from the solid wall for supersonic flow 
over a flat plate. 
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Figure 5.10. Surface pressure distribution for the viscous internal flow over a 10-degree 
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b. CONCLUSIONS 


The second-order PNS algorithm due to Schiff and Steger (Ref, 1) has 
been coded and run successfully for three model problems. Good results are 
obtained for the flow over a flat plate with an impinging shock wave. Solu- 
tions of the flows over a 10-degree wedge indicate the difference in the 
peak pressure locations. This difference is believed to be due to the 
splitting of the shock waves which, in turn, is caused by oscillations of 
the solutions. Investigations are being made to improve the peak pressure 
1 ocations. 

Also, the formulation of a fairly general second-order marching algo- 
rithm is under current investigation. This algorithm will not be restricted 
to only a three-point-backward differencing scheme of Ref. 1. The results 
obtained from this general algorithm will be compared with other PNS codes 
as well as with the full Navier-Stokes solutions. 
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